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I.  INTRODUCTION 


This  thesis  proposes  an  improvement  to  present  near-shore  wave  predietion 
models.  Aeeurate  predietion  of  wave  shoaling  transformation  is  erueial  for  the  sueeessful 
planning  of  amphibious  operations,  providing  a  strong  military  motivation  for  this 
researeh.  Additionally,  the  aetion  of  waves  and  wave-driven  eurrents  on  a  beaeh  drives 
sediment  transport  and  beaeh  erosion,  stimulating  eivilian  interest  in  this  field  as  well. 

As  oeean  surfaee  waves  shoal  from  deep  to  shallow  water,  amplitudes  inerease, 
wavelengths  deerease,  and  propagation  direetions  refraet  toward  normal  ineidenee  to  the 
beaeh  as  a  result  of  linear  proeesses.  Additionally,  non-linear  interaetions  ean  eause 
signifieant  energy  transfers  to  wave  eomponents  with  both  higher  and  lower  frequeneies 
than  the  original  waveform.  In  deep  water  (  kA  »  1 ,  where  k,  is  the  wave  number  and  h 
is  the  water  depth)  and  intermediate  depths  {nh  =  0(l)),  nonlinear  effeets  on  waves  are 

usually  evaluated  with  finite  depth  theory  on  the  basis  of  Stokes  perturbation  expansion 
for  small  wave  slopes  (aK«l,  where  a  is  the  wave  amplitude)  (e.g.,  Phillips,  1960; 
Hasselmann,  1962).  Nonlinear  interaetion  between  two  primary  wave  eomponents 
exeites  seeondary  waves  with  frequeneies  and  wave  numbers  that  are  the  sum  (or 
differenee)  of  those  of  the  primary  waves.  In  deep  and  intermediate  water  depths,  these 
"triad"  interaetions  are  nonresonant,  and  the  resulting  seeondary  wave  amplitudes  are 
small  eompared  to  the  primary  wave  amplitudes  {Phillips,  1960;  Herbers  et  al,  1992, 
1994).  Similarly,  tertiary  waves  may  be  exeited  by  non-linear  interaetion  between  three 
primary  wave  eomponents,  eonstituting  "quartet"  interaetions.  For  eertain  eombinations 


of  wave  numbers  quartet  interaetions  are  resonant,  eausing  a  gradual  exehange  of  energy 
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between  the  wave  eomponents  that  is  important  to  the  evolution  of  wind-wave  speetra 
over  large  distanees  I960;  Hasselmann,  1962;  Hasselmann  et  ai,  1973). 

Unlike  the  deep  and  intermediate  water  regimes,  in  shallow  water  (  k/?  «  1 )  the 
seeondary  response  is  strongly  amplified  due  to  near-resonanee.  The  eonventional 
approaeh  to  deseribe  these  interaetions  is  the  use  of  Boussinesq-type  equations,  whieh  are 
based  on  the  assumption  that  measures  of  nonlinearity,  e  {=a/h),  and  dispersion 

cr^(=  (nhY)  are  both  small  and  of  the  same  order.  This  approaeh  reduees  the  three- 
dimensional  fluid  motion  problem  into  a  two-dimensional  one  by  introdueing  a 
polynomial  approximation  of  the  vertieal  strueture  of  the  flow  field,  while  aeeounting  for 
weak  non-hydrostatie  effeets  due  to  the  vertieal  fluid  aeeeleration  (see  Madsen  & 
Schaffer,  1998,  for  a  eomprehensive  review).  Boussinesq’s  original  work  (1872)  was 
restrieted  to  a  horizontal  bottom;  sinee  then  that  work  has  been  extended  to  an  uneven 
bottom  in  one  dimension  {Mei  &  LeMehaute,  1966),  and  later  two  dimensions 
{Peregrine,  1967). 

Berbers  and  Burton  (1997,  referred  to  as  HB  in  the  following)  formulate  shoaling 

evolution  equations  for  the  speetral  amplitudes  and  phases  of  a  direetionally  spread  wave 

field  propagating  over  a  shallow,  gently  sloping  beaeh  with  straight  and  parallel  depth 

eontours  that  are  aeeurate  to  the  lowest  order  ( e  )  in  nonlinearity  and  dispersion  (see  also 

Freilich  &  Guza,  1984;  and  Elgar  &  Guza,  1985).  Limitations  of  HB  inelude:  1)  negleet 

of  alongshore  depth  variation,  2)  restrietion  to  gently  sloping  beaehes,  3)  restrietion 

to  small  oblique  wave  angles,  4)  a  low  order  approximation  of  dispersion  effeets  whieh 

limits  the  model  to  very  shallow  water  (nominally  depths  less  than  10m),  and  5)  negleet 

of  nonlinear  interaetions  between  more  than  three  wave  eomponents.  In  this  paper,  we 
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propose  improvements  that  address  the  limitations  mentioned  above  by  earrying  out  the 
HB  derivation  to  one  higher  order  ( )  in  nonlinearity,  dispersion,  and  bottom  slope,  and 
allowing  for  weak  alongshore  depth  variations.  This  higher  order  model  still  assumes  a 
gently  sloping  beach  and  small  oblique  angles,  but  improves  the  accuracy  of  the  model 
for  finite  bottom  slopes  and  angles.  Dispersion,  directional  spreading,  and  strength  of 
triad  interactions  are  all  represented  to  one  higher  order,  as  well  as  the  effects  of  the 
cross-shore  bottom  slope.  New  effects  introduced  at  the  lowest  order  are  alongshore 
depth  variation  and  nonlinear  quartet  interactions. 

An  overview  of  the  technique  used  to  obtain  this  result  is  as  follows;  The  first 
order  horizontal  flow  and  pressure  is  shown  to  be  depth  independent.  The  continuity 
equation  is  then  depth-integrated  to  obtain  a  lowest  order  expression  for  the  vertical  flow. 
This  vertical  flow  solution  is  substituted  into  the  irrotationality  condition  and  the  vertical 
momentum  equation,  which  are  then  depth-integrated  to  obtain  the  next  higher  order 
expressions  for  the  horizontal  flow  and  pressure,  respectively.  This  process  is  repeated 
until  third  order  expressions  for  the  horizontal  flow  and  pressure  are  obtained.  Following 
Freilich  and  Guza  (1984)  and  HB,  an  assumed  first  order  spectral  wave  field  is 
substituted  in  the  governing  equation  for  the  depth  averaged  velocity  potential,  and  the 
solvability  condition  yields  a  coupled  set  of  first  order  equations  for  the  evolution  of  the 
spectral  amplitudes. 

A  significant  challenge  to  arriving  at  such  higher  order  solutions  is  that  the 
complexity  of  the  polynomials  at  each  level  of  the  derivation  grows  significantly  as 
higher  orders  of  non-linearity  and  dispersion  are  considered.  Madsen  &  Schaffer  (1998), 
Agnon  et  al.  (1999),  and  Gobbi  et  al.  (2000)  derive  a  variety  of  Boussinesq-type 
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equations  to  high  orders  of  accuracy,  from  which  this  complexity  is  quite  evident. 
Although  equations  derived  in  these  earlier  studies  are  more  general  than  those  that 
follow  in  this  thesis,  they  are  not  cast  in  a  spectral  form  and  require  prohibitively 
expensive  full  time-space  domain  solution  techniques. 

This  thesis  is  organized  as  follows.  In  Chapter  II  the  formulation  of  governing 
fluid  flow  equations  are  expanded  to  third  order  and  integrated  over  the  water  column.  A 
coupled  set  of  amplitude  evolution  equations  for  a  spectrum  of  waves,  the  main  result  of 
this  paper,  is  derived  in  Chapter  III.  In  Chapter  IV,  linear  terms  in  these  equations  are 
verified  by  analytical  methods,  and  example  computations  for  a  plane  beach  quantify 
some  of  the  improvements  of  this  model  over  the  lower  order  HB  model.  Finally, 
Chapter  V  summarizes  these  results,  and  identifies  opportunities  for  further  research. 
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II.  THIRD  ORDER  GOVERNING  EQUATIONS 


A.  ASSUMPTIONS 

The  equations  of  motion  and  the  eontinuity  equation  may  be  written  in 
nondimensionalized  form  (neglecting  viscous  effects): 


— +  (u*V)u  +  w  — +  Vi?  =  0, 
ot  oz 

dw  /  dw  dp 

- |-(u*V)w  + w - h  — +  1  =  0, 

dt  dz  dz 


V*uH - 

dz 


0, 


(1) 


where  t  is  time,  V  denotes  the  two-dimensional  gradient  operator  {d ! dx,d ! dy),\\  is  the 
vector  (m,v)  of  the  horizontal  {x,y)  velocity  components,  w  is  the  vertical  (z)  velocity 
component,  and  p  is  pressure.  The  x-axis  points  onshore,  y  points  alongshore,  and  z 
points  upward  with  z  =  0  corresponding  to  the  mean  sea  surface.  The  surface  and 
bottom  are  defined  by  r](x,y,t)  and  z  =  —h{x,y)  ,  as  shown  in  Figure  2.1. 


z. 


Figure  2.1.  Definitional  sketch  of  variables  and  coordinate  frame.  Directionally  spread 
waves  propagate  over  a  beach  with  gentle  cross-shore  and  alongshore  depth  gradients. 
The  variable  6  denotes  the  wave  propagation  direction,  k  is  the  wavenumber  vector,  and 

/zq  is  a  representative  water  depth. 
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All  variables  in  equation  |(1)|  are  normalized  using  gravity  g,  density  of  seawater 
p ,  and  a  representative  water  depth  .  Irrotational  flow  is  assumed,  so  that 


du  _  dw  dv  _  dw  du  _  dv 

dz  dx  ’  dz  dy  ’  dy  dx 

The  surfaee  and  bottom  boundary  eonditions  are  given  by 

dri  _ 

w  =  — +  u*Vr/,  z  =  T], 

at 

T*  =  0,  z  =  r}, 

w  =  — u*V/2,  z  =  —h. 


(2) 


(3) 


The  variables  p,  u,  v,  w,  and  p  are  expanded  in  terms  of  the  non-linearity 
parameter  e,  defined  to  be  the  ratio  of  wave  amplitude,  a,  to  water  depth,  h.  In  shallow 
water,  v  and  w  are  both  smaller  than  u  (where  k  is  the  wavenumber),  and  thus  the 

perturbation  expansions  of  the  sealed  variables  are  given  by 

U  EU^  E  II2  E  W3  “h 

V  =  a(ev^  +  + ...),  (4) 

w  =  <t(£Wj  +  e^W2  + 
p  =  -z  +  epi  +  P2  +  E  P2  + 


where  a  =  nh  is  a  measure  of  dispersion.  The  independent  variables  are  then  sealed 


d  _  d  d  _  2  d  d  _  d  d  _  d 

dx  dx'  ’  dy  dy'  ’  dz  dz'  ’  dt  dt 


so  that  the  order  of  the  terms  appears  explicitly  when  equations  |(4)Jand|(^are  substituted 
into  the  governing  equations.  In  the  classic  Boussinesq  approximation,  dispersion  E 


and  nonlinearity  e  are  assumed  to  be  of  the  same  order.  The  water  depth  is  taken  to  be  a 
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function  of  the  slow  variables  (x.y)  sueh  that  -^^=3^  and  =  where 

dx'  ^  dx  dy'  dy 

(3  =  Vh  /  nh  is  assumed  to  be  0[e).  In  this  approximation,  alongshore  depth  variations 
are  smaller  than  eross-shore  variations. 

B,  VERTICAL  STRUCTURE  OF  THE  FLOW  AND  PRESSURE  FIELD 

In  this  seetion  the  irrotationality  eonditions,  vertieal  momentum  equation,  and 
eontinuity  equation  are  integrated  over  the  water  eolumn  at  sueeessive  orders  of  e  to 
determine  the  vertieal  strueture  of  the  flow  and  pressure  field  aecurate  to  First, 

we  expand  the  irrotationality  eonditions  1^2)  [and  vertical  momentum  equation  (lb)  using 


(4b-e)  andl(5)jto  0{e^y. 

du, 


dz 

cIVj 


j  +  e 


dz' 


+  s 


dU2 

dwA 

[dz' 

dx'j 

[dV2 

dw,' 

dz' 

dy'] 

+  r 


du^ 

dw^] 

[dz' 

dx'] 

dv^ 

dw^] 

dz' 

dy'\ 

=  0(6^), 


(6) 


(V) 


du^ 

> 

1  c 

dU2 

dv2 

du^  dv^ 

w 

dx' 

/ 

c 

w 

dx' 

/ 

c 

dy'  dx'  ^ 

Oi£^), 


(8) 


^  I 

dp  2  dw,' 

dz' 

r  c  1 

[dz'  dt'] 

+  £ 


dp  2 

,dw2 

r 

h 

[dz' 

^  dt' 

w, 


dz' 


=  0(e^). 


(9) 


Note  that  ui,  vi,  and  pi  are  independent  of  z'  to  0[e).  Substitution  of  the  dynamie 
surfaee  boundary  condition  (3b)  in  (4e)  yields 

Pi=Vv  (10) 
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Next  we  expand  the  eontinuity  equation  (le)  and  the  kinematie  bottom  boundary 


eondition  (3o)  and  eolleet  terms  of  like  order, 


du^ 


0, 


Wj  =  0  at  =  —h. 


(11) 


du.  dv,  dw-, 

— - - ^  =  0 

dx'  dy'  dz' 


dh  ! 

w,  =  —u,  at  z  =  — , 
ox 


At  the  lowest  order,  equation  |(1 1)  [integrates  to 


(12) 


Wi 


=  -(z'  +  /z) 


(13) 


Substitution  of  wi  and  pi  into  [7)J  and  [9)  [and  integration  over  the  water  eolumn 
yields  the  vertieal  distribution  of  the  seeond-order  veloeity  and  pressure  field 


^2  ^2 


Iz'^ 


+  hz' 


d^u. 


dx 


n  ’ 


(14) 


V2  =^2 


‘  z'^ 

- \-hz' 

2 


dx'dy'  ’ 


(15) 


P2=  Pi  - 


‘  z'^ 

- \-hz' 

2 


dx'dt' ' 


(16) 


where  U2,  V2,  and  P2  are  integration  eonstants. 


The  eonstant  P2  follows  from  substituting  p\  and  p2  in  |(4)|  and  evaluating  the 
dynamie  surfaee  boundary  eondition  (3b)  to 


Pi  —  Vi- 


(17) 
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The  remaining  eonstants  U2  and  V2  can  be  expressed  in  terms  of  the  second-order  depth- 


averaged  flow  variables  Wj  and  V2  tequations  It73tlandlt76tlin  Appendix  A).  With  these 

results,  the  second  order  wave  motion  can  be  written  as 


,  /  h 

- \-hz  — 

2  3 


^  ^  d^u, 


dx 


n  ’ 


(18) 


- \-  hz  -\ - — .  . , 

2  3  ]dx'dy' 


(19) 


Pi  =V2 


+  hz' 


^  d^u. 


/ci./  • 


dx  dt 


(20) 


Having  determined  the  second  order  pressure  and  horizontal  flow  variables,  we 


may  proceed  as  in  the  case  of  wi  and  integrate  |(12)|to  determine  W2: 


'z'^ 

hz'^ 

1 

h^z'] 

1 

d^’Uy 

-(7^  -\^h\ 

du^ 

> 

j 

6  ’ 

r 

2 

r 

3 

/ 

dx'^ 

—  1  z  -r  n\ 

dx' 

"a/I 

dh 

—  u, 

dx 


(21) 


Next  we  evaluate  1/3  by  substituting  wi,  U2,  V2,  and  W2  into  equation [(^  and  integrating  the 
terms  with  respect  to  z': 


■■u,- 


(  n 

^  -L  hP 

d^U2 

7^ 

CM 

CO 

J 

+ 

(  /4 

Z 

hz"" 

,2  I2\ 

h  z 

- r  fiz. 

2 

dx'^ 

dx'dy' 

r 

6 

r 

6 

d‘^u. 


dx 


/4 


dh  du. 

-2z^ — 7, 

dx  dx 


(22) 


Expressing  the  integration  constant  U3  in  terms  of  the  depth-averaged  third  order  flow 
using  equation  |(74)  [(Appendix  A),  we  obtain 
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M,  =  M,  + 


hz 


/3 


2_/2 


h^Z 


24 


45 


c)x 


/4 


- \-hz  -\ - 

2  3 


d^U2 

o2  ) 

1  d  Vi 

dx'^ 

V 

dx'dy' 

-(2z'  + a) 


,  I  ,  \  hrj^ 

dx  dx'  3  5x 


/2  • 


(23) 


Similarly  for  vs, 


V3=F3- 


‘  z'^ 

- \-hz' 

2 


^  1  a2 


c)  M,  d  V. 
2  ^  1 


2_.  1  r_/4 


dx  dy  dy 


n 


+ 


hz'^  z^  1 


24 


dx'^dy' 


—z 


dh  du^ 

dh  du, 

1  1 

dy  dx' 

dx  dy' 

(24) 


and  substituting  for  F3  using  equation |(7^ (Appendix  A)  yields 


+ 

1 

II 

f  z'"  hz'^ 

/2V'  A"! 

d'^u^ 

1 

24  '  6 

1 

6  45 

/ 

dx"’dy' 

^  2 

h^\{  d^U2  d^Vi 


dx  dy  dy 


n 


,h 

z  +  — 


duy  dh  duy  dh 


dy'  dx  dx  dy] 


rjji  d^Ui 


/a../  • 


3  dx  dy 


(25) 


Finally,  pj,  is  determined  by  substituting  p\,  p2,  wi,  and  rv2  into  the  vertical 
momentum  equation  (lb)  and  depth-integrating  the  0(^e^]  terms: 


Pi=Pi  + 


( z'^ 

- \-hz' 

2 


d^u. 


d'^Vy  d^Uy 


dx'dt'  dy'dt'  ‘  dx''^ 


du^ 


dx' 


^  z'""  hz 

- 1 - 

24  6 


/3 


d^Uy  ,  dh  duy 


(26) 


dx'^dt'  dx  dt' 


From  equation  (3b)  it  follows  that 
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(27) 


(^1  +  £7?2  +  £^^3  +  0{e^  ))  +  sVi 

j'  o2  'j 

+e^  %  +  erj^h  +  0(£^)  +£^  (Pj  +  0(6:))  + 0(£'^)  =  0. 

Hence, 

O^u 

P  -r^  r],h-=^  (28) 

'  ''  "  dx'dt' 

and 


A 


dx'dt'  ^ 


ay'a^' 


+  Mj 


z''^  hz'^  z'^  d'^u,  !  dh  du,  ,  d^u, 

24  6  6  j  dx  ^dt  dx  dt  dx  dt 


(29) 


In  summary,  the  flow  and  pressure  fields  are  described  by  equations  |(10)  [and 
(13)-  (29)j  Equations  1(10)  land  1(1 3)H(20)]are  identical  to  those  in  HB.  The  higher  order 
equations,  |(21)||(29)J  improve  the  HB  results  by  expressing  the  vertical  flow  structure  as  a 
fourth  order  polynomial  in  depth,  and  by  including  nonlinear  and  bottom  slope  terms  that 
are  neglected  at  lower  orders. 

C.  REDUCTION  OF  THE  GOVERNING  EQUATIONS 

Here  the  velocity  and  pressure  distributions  derived  in  the  previous  section  are 
substituted  into  the  governing  equations  to  obtain  depth-integrated  horizontal  momentum 
and  continuity  equations.  First,  the  perturbation  expansions  |(4)  |and  scaling  relations  (5) 
are  substituted  into  the  horizontal  momentum  equation  (la).  Collecting  terms  of  like 
order. 
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du.  dp, 
dt'  dx' 


+  s 


^  duj  du.  dp  ^ 

_ ^-1-77  _ Lj _ £_ 

dt'  ^ ' ' 


dx'  dx' 


du^ 

,  a(MiM2) 

du, 

1-11  i  1 

dt' 

'■  dx'  ^ 

du^ 

J 

dz' 

L 

0{^) 


and 


.  dPi 
dt'  dy' 


+  e 


dv. 


dt 


T  +  ^i 


avj 

,dP2 

dx'] 

+  £ 


dv. 


dv. 


dv 


dt'  '  dx' 


dx 


1 

'y  +  y  ' 


dv 


dy'  '  dz'  '  dy' 


2  ^  dp^ 


=  0{e\ 


Substituting  wi,  U2,  V2,  pi,  W3,  V3  andp^  (equations  |(  13)|  |(  1 8)H(20)i  |(23)1  |(25)i  |(29)|)  yields 
after  some  algebraic  manipulation 


duy 

,  A 

1—  c 

du2 

,9% 

du, 

1  7/  ^ 

A"  a\ 

dt' 

ax' 

r  c 

dt' 

dx' 

^  'dx' 

3  dx'^dt' 

+£ 


du^ 

,  A 

A" 

d^u. 

2  1 

a"vj 

h*  d^u, 

dt' 

ax' 

"  3 

dx'^dt' 

dx'dy' dt' 

45  dx'^dt' 

+ 


h  d 

1>W 

2 1 


2..  t 


Vi 


d  Uy 


dx 


n 


d  /  — 


dx 


+  V, 


du^ 

dy' 


h 

H - 

3 


.  du,  d  u, 

2 - 7 - TT  —  U, 

dx'  dx'^  ' 


1  ,  d 


dx 


n 


-h 


dx' 


Vi 


d^u^ 


dx'dt' 


h 


dh  d^u^ 
dx  dx'dt' 


0{e% 


(30) 


A  1  - 

av2 

at'  ay' 

at' 

u 

5vj  d^u^ 


+£ 


av3  dp2 

a^Mj  a"vj 

d^u, 

h 

dh  d^u^ 

dh  d^u, 

1  1 

dt'  dy' 

"  3 

dx'dy' dt'  ^  dy'^dt'  ^ 

45  dx'dy  dt' 

’2 

dy  dx'dt' 

dx  dy'dt' 

^3  at' 

+  h 


Vi 


d^u. 


dx'dy' 


+  u. 


dv. 


dx' 


■  dv,  dv, 

— ^  +  Vj — ^ 

dx'  dy' 


du^  d^u^ 

h^  a 

dx'  dx'dy' 

3  ay'i 

2..  t 


a  Mj 


'ax'" 


-h 


a 

ay' 


Vi 


d^u, 


dx'dt' 


0{£^). 


(31) 
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Finally,  integrating  the  continuity  equation  (Ic)  over  the  entire  water  column, 
applying  Leibniz  Theorem,  and  the  kinematic  surface  and  bottom  boundary  conditions 
(3a,c), 


ry  ry  H 

or]  o  r  r  or, 

- 1 - /  uaz  H - 7  /  vaz  = 

dt  dx  ^  ^ 


d_ 


(32) 


Substituting  the  expansions  (4a-c)  and |(5)| yields 


d 


-f{srj^  +s^r]2  +0(£"^))  +  cr — j  J {eu^  +  ^^^3  -\-0{e^)^dz' 


dt'^  "  ^  dx 


—  h 


o  'I 

+cr^ — -  J  a(^ev^  +  e^V2 +0(s^  )'jdz' =  0. 


(33) 


Now  evaluate  the  integrals  for  each  of  the  flow  variables  using  equations  (4a),  |(18)|  |(19) 


K23)  andpS) 


J  Uydz' =  {h  +  eri^-\-e^r]^u^-\-0{e^). 


(34) 


—h 


'I 

I 

—h 


U2dz' 


[h  +  erj, )  U2  -  er],h^  +  0(6^ ), 


(35) 


J  u^dz'  =  hu^  +  ^  ^2  +0(5),  (36) 

3  dx 


V 

J  v^dz'  =  {h-\-  £?7i )  Vj  +  0{e^ ), 

—h 


(37) 


V  _ 

j  Vj  =  hv2  +  0{e). 

-h 


(38) 
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Combining  and  ordering  the  terms, 


%  +  *^  +  e 

dt  dx 


+£ 


drii 


dt' 


f  A 


dU2 

dv, 

d 

1 

dx' 

dy', 

dx' 

dVi 

h 

du^ 

dv2 

^  d[viU2) 

dt' 

r 

dx' 

dy'\ 

^  dx' 

d  ,  X  — dh  dh 


dx' 

0(.=  ). 


dx 


(39) 


The  system  of  equations  |(30)J  |(31)J  and  |(39)|  eonstitute  the  full  set  of  governing 
equations  to  O  .  This  system  is  now  eross-differentiated  to  eliminate  773 ,  whieh  gives 


d 

d 

de 

dx 

du^ 


dx 


+£ 


d^U2  d 

/ 

h 

dU2 

dv; 

d 

h  d^  [hu^) 

df  dx 

n 

dx 

dy 

dx 

3  dx^df 

A 

dx 


dh 

[d~x 


+ 


d 

dt 


du. 


dx 


dx 


+£ 


d^u^ 

d 

/ 

h 

du^ 

1 

1 

CO 

J 

\  0 
h^ 

d‘^U2 

df 

dx 

n 

dx 

^  6y] 

3 

dx^de 

d\ 


dxdydt 


d%, 


dx 

d 

^  dt 


dh—  dh 


— —u^  +  ■ 
dx 


dy 


-h 


d^ 


du2  —  du^ 


h^  d 


dx 

( 


dx 


d 

^  dt 


dxdt 
( 


Vi 


d^u^ 

h  d^ 

1 

dxdt 

/ 

zdf 

45  dx^de 
2. 


Vi 


d 


dx^ 


du^ 

h^  d^ 

duy 

dy , 

1 

2  dxdt 

,  dx  , 

3  dxdt 


u, 


d^u. 


dx^ 


d^  /  - 


_  0 

=0(£^) 


dx^ 


dxdy 


(40) 


and 
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where  the  primes  are  dropped  from  here  on  to  simplify  the  notation. 

D,  GOVERNING  EQUATIONS  IN  VELOCITY  POTENTIAL  FORM 

Finally,  expressing  the  horizontal  flow  variables  as  gradients  of  veloeity  potential 

functions  0,  (see  Appendix  B)  allows  equations  |(40)|  and  1(41)1  to  be  converted  to  the 
gradient  form 
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Since  equation  (42)  must  hold  in  the  entire  horizontal  plane,  it  follows  that  the 
expression  inside  the  gradient  operator  (disearding  a  uniform  time  dependenee)  must 
vanish.  Therefore,  the  final  form  of  the  governing  equation  is 
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de 


+£ 


dx^ 


dt  [  dx^ 
dh  c)0j  1  d 


dy^ 

d(j), 


aVi 


3  dx  of 


dx  dx  2  dt 


^2 


dx 


A 

c)x 


501 


c)x 


+£ 


5^3 

-  h 

d^(j), 

5^3  1 

h^  d^ 

\d^cj>2 

,  d\  ] 

h^  aVi 

df 

n 

dx^ 

dy^  \ 

3  df 

dx^ 

d/\ 

45  dx^df 

dh  d(j)2  dh  dcf)^  .  dh  d^cf)^  .  c) 


h^^^^-h  — 
dx  dx  dy  dy  dx  dxdf  dt 


Vi 


dx  dt 


d 

^  dt 
dx 


dcj),  d^2 


dx  dx 

d^2 


1  d 
^ - 

dcj), 

'  h^  d 

H - 

/  ■>  \ 
d^A 

2  0 
h^  d 

d(f>i 

2dt 

dy 

2  dt 

dx^ 

\  2 

3  dt 

dx  dx^ 

Vi 


dx 


d 

L  ^<^01 

+  n 

OX 

7 

dy 

/  dy] 

+  O(e^)  =  0. 


(43) 
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III.  THIRD  ORDER  AMPLITUDE  EVOLUTION  EQUATION 


A,  ASSUMED  SPECTRAL  WAVEFORM 

In  the  previous  chapter,  we  derived  a  governing  equation  |(43)|  for  the  depth- 

averaged  velocity  potential  function  to  0(£^^.  We  now  derive  spectral  amplitude 
evolution  equations  by  assuming  the  lowest  order  wave  field  to  be  a  linear 

superposition  of  nearly  plane,  shoreward-propagating  waves 

(^’  y)  (^’  ’  0]  ’  (44) 

p  ^ 


\-  I  ^\ 
=  - exp  {x,y,t) 


where  =  p^uj  and  =  ^A/  are  the  scaled  frequency  and  initial  alongshore 

wavenumber,  with  Ao;  and  A/  the  separation  of  adjacent  bands  in  the  Fourier 
representation.  The  phase  function  ^  contains  the  fast  variations  of  the  lowest-order 

wave  field. 


X 


The  complex  amplitude  is  a  slow  function  of  x  andy 


+  o (e" ) , 


(46) 


(47) 


where  R\p,q,  Rip.q,  and  S\p_q  account  for  amplitude  and  phase  changes  resulting  from  linear 

(e.g.,  dispersion,  shoaling,  refraction)  and  nonlinear  (triad  and  quartet  interactions) 

19 


processes,  and  ehanges  in  assoeiated  with  alongshore  depth  variations  are  absorbed  in 
the  slow  alongshore  variations  of  It  will  be  shown  that  dA^  ^  jdy  terms  appear  at 
suffieiently  high  order  that  they  do  not  require  explieit  expression  beyond  0(e) . 

B,  0(1)  PROBLEM 

Define  the  linear  operator  X  as 


(48) 


It  is  easily  verified  that  is  a  solution  to 

£[4>,)  =  0{e).  (49) 

To  determine  solutions  to  higher  order  problems,  the  residual  right  hand  side  of|(49)|must 
be  evaluated.  To  0{e^^,  is  given  hy  equation  1(91  )lin  Appendix  C 

C.  0{e)  PROBLEM 

At  O  (e) ,  applying  the  results  of  Appendix  C, 


2R  A  +i 

j  — 

p  ? 

3 

r 

2h 


p—m^q—n 


(50) 


exp  +  0{e). 


Resonant  growth  of  02  is  prevented  by  the  solubility  eonstraint  that  the  right-hand  side 


foreing  terms  in  (50)  do  not  eontain  any  free  wave  solutions  of  the  general  form 
.  Sinee  all  these  terms  obey  to  0[e)  the  shallow  water  dispersion  relation. 
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it  follows  that  the  right  hand  side  of|(50)|must  vanish.  Solving  for  R\p_q  yields  the  lowest 


order  term  in  the  order  amplitude  evolution  equation: 

R\p,q  =  ^lp,q^p,q  ^  Q\ p'Y^'Y^  A ,A p-r q 

r  s 

where 


-Zl^  ■ 
~  4hdx^^ 


(51) 


6 

2^J 

3uj 


This  result  was  derived  in  HB. 

D,  0[e^)  PROBLEM 

To  determine  the  remaining  higher  order  amplitude  evolution  terms  R2p,q  and  Sip_q, 
we  proceed  in  the  same  fashion  as  for  Rip,q.  Since  xi^(p2j  =  0[£),  represents  (to 

lowest  order)  free  shallow  water  wave  motion  that  can  be  absorbed  in  (/>,.,  and  thus  (^>2 
can  be  set  equal  to  zero.  Note,  however,  that  rj^  (evaluated  in  Appendix  D)  contains 
deviations  from  free  shallow  water  wave  motion  that  contribute  to  the  forcing  of  . 


Tedious  but  straightforward  evaluation  of  all  the  forcing  terms  in  equation  1(43) 

(see  Appendix  C)  yields 
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Z((/>3l  = 


EE 

P  q 
\{ 


2hl. 


ih  dR 


^p,‘i 


2p,?  Ip,? 

S  dx 


2h^^^ujl  i  dh 


u)p  dx 


R 


i/',? 


+ 


dh 

_  h  dh 

6 

dx 

^P  dy^ 

+  i 


^  -h^Lvl  h^LOjl^ 


45 


A 


P,<i 


(52) 


4,.v  dh 

,  24..  1 

1  7* 

5u)pUJ^  to 

A 

2h^d 

h^d 

/ 

1  1 

“r  0 

[3  4  J 

21.  ,v 

+EE 

r  s 

+EEEE 

r  s  m  n 

Finally,  substituting  (equation  |(51)||  and 


A 


p—r,q—s 


-iUJp 

AAA 

m,n  r,s  p—r—m,q—s—n 

dR 


-1  d^h 


dx 

+EE 

r  s 

+EEEE 

r  s  m  n 


Ah  16 


dx 


+  h 


-^+ 

.2  r2 

f 

36  ^ 

6 

H] 

34 

( 

1 

Ah 

2 

\  ^ 

_ £_ _ J_ 

3  ^pj 
9<^p^p-r 


+ 


(  3 

/n 

r 

3 

+  i- 


21u}„  dh 


I6h^^^  dx 

Am,nA,s\-r-m,q-s-n  +  ^  (e) 


A^ 


A  A 

'^r,s'^p—r,q—s 


in  equation  (52)  and  setting  the  forcing  terms  on  the  right  hand  side  equal  to  zero  to 


remove  all  secular  terms,  the  following  relation  for  the  higher  order  amplitude  gradients 


R2p,q  and  Sip^q  is  obtained: 


h^'^l. 


^2p,q  +  -  ^2p,q^p,q  +  ^'^Q2p,q,r,sA,s^p-r, 

r  s 

+  EC2,,,EEE4.,4,^ 


(53) 


p—r—m,q—s—n  ’ 


s  m  n 


where 


22 


L 


'2p,q 


~4 

-hi 


Cl 

dh 

h 

dh 

dx 

2h''\ 

dy 

1  C  1 

1 

1 

idh] 

'  1  dA 

360  12 

.dx, 

d7 

Q: 


2p,q,r,s 


-15  dh 
32h^  dx 


+  i 


^2p,r  ~  ^ 


+ 


-5Up  +  9(jO^ 

16^^ 


3 

f  .3  72; 

3 

/M 

16/2*/' 

3 

8*''- 

3 

1  p 

r  / 

+ ^^pA 

^  - 

12/2'/" 


Ih^'^uol 


E,  EVOLUTION  EQUATION 

Combining |(47)[  [(STjj  and  (53),  our  final  result  in  amplitude  form  is 


dA„  „  dA 


dx 


+ 


<^p  dy 


hlp,g  +  Lip,^  )  +  ^^{Qlp  +  Q2p,q,r,s  )  A, A, 


'p—r,q—s 


+Ec2,„EEE4,.4,,-4, 


'p—r—m,q—s—n 


+  0(s--). 


(54) 


s  m  n 


The  real  part  of  the  linear  term  on  the  right  hand  side  is  the  amplitude  growth  of  the  given 
mode  due  to  shoaling  (i.e.,  ehanges  in  the  group  speed);  the  imaginary  part  represents 
slow  phase  ehanges  owing  to  dispersion,  direetional  spreading,  and  bottom  slope  effects. 
The  double  and  quadruple  sums  contain  the  slow  amplitude  and  phase  changes  of  a  mode 
resulting  from  nonlinear  interactions  of  all  possible  triads  and  quartets  (respectively)  in 
which  the  mode  participates. 

The  new  terms  (relative  to  HB;  L2p,q,  Qip.q.r.s,  and  Cip.r)  introduce  improved 
representations  of  dispersion,  directional  spreading,  triad  interactions,  and  the  effects  of 
the  cross-shore  bottom  slope  (all  to  one  higher  order  of  accuracy),  as  well  as  the  lowest 
order  effects  of  alongshore  depth  variations  and  nonlinear  quartet  interactions  not 
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represented  in  HB.  The  verification  of  many  of  these  new  terms,  and  a  quantitative 
investigation  of  the  improvements  gained  by  them,  is  given  in  the  next  chapter. 
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IV.  VERIFICATION  OF  LINEARIZED  MODEL  FOR  AN 
ALONGSHORE  UNIFORM  BEACH 


A,  VELOCITY  PROFILES 

From  linear  finite  depth  theory,  the  horizontal  velocity  profile  (in  the  x- 
direction)  can  be  expressed  in  normalized  form  as 


u 

u 


cosh(«;(z  +  /2))  «:Acosh(«;(z  +  /2)) 


cosh(K(z +  /z))  Jz 


smh(reA) 


(55) 


Expanding  the  hyperbolic  functions  for  small  nh  and  using  the  geometric  expansion  for 
the  quotient  in  equation  1(55)  yields 


nh 


u 
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[z  +  hf  + 
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■O^nhy 
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Collecting  and  combining  like  terms, 


u 

u 


1  +  k" 

(^2  7  2) 

z  ,  n 

- \-hz-\ - 

[2  3j 

^  z‘^  hz^  h^z 


3  ,.2_2 
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45 


+  0[Khf 


(56) 


The  z -polynomials  at  each  order  in  equation  (56)  are  in  exact  agreement  with  the 


polynomial  coefficients  at  the  corresponding  orders  for  U2,  V2,  M3,  and  V3  (equations  [18) 


iml  t23}J  and  to).  Figure  4.1  compares  vertical  structure  of  the  first  (mj  ) ,  second 
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(mj  +  £M2  )  5  third  (u^  +  eu^  +  )  order  approximations  of  u  with  the  exact  linear 


finite  depth  solution  |(55)|for  a  wave  with  period  10s  in  20m  water  depth.  The  first  order 


approximation  is  the  familiar  depth-independent  shallow  water  solution  for  the  horizontal 
flow.  The  second  order  approximation  better  represents  the  depth  dependence  of  the 
horizontal  flow,  while  the  third  order  approximation  virtually  replicates  the  exact  finite 
depth  solution. 
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Figure  4.1.  Vertical  profiles  of  normalized  horizontal  current .  Note  that  the  first  order 
estimate  is  depth  independent.  Higher  order  estimates  converge  toward  the  vertical 
structure  predicted  by  finite  depth  theory. 


B,  DISPERSION 

In  linearized  form  ,  the  amplitude  evolution  equation  for  an  alongshore 
uniform  beach  {dhjdy  =  0,  dA/dy  =  O)  is  given  by 
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(1-1 

1 

r2  Y 

—  +  - 

4 - V 

11 4* 

4 

[  ^l\] 

'  P  '  ^  I  /^3/2 

6 


11^: ,  c  ] ,  1  [dh^ _ L 

360  12  86jM  32/2'/'cu  laxj  8/2‘^'w 


The  imaginary  part  of  the  right  hand  side  represents  slow  phase  ehanges  as  the 
wave  eomponent  evolves  in  the  eross-shore  direetion,  that  is,  eorreetions  to  the  eross- 
shore  wavenumber,  k.  From  equation  '57)  it  follows  that  k  for  a  wave  with  radial 

frequeney  u  and  alongshore  wavenumber  /  in  uniform  depth  is  given  to  (9  (re  h)^  by 


nh^_h^ 

~  ^  6  ^  360  2u  ^  12 


Now  eonsider  the  exaet  k  in  finite  depth.  The  alongshore  and  eross-shore 
wavenumbers  satisfy  the  relationship 

k^+l^=K\  (59) 

Using  the  finite  depth  dispersion  relation,  the  eross-shore  wavenumber  ean  be  expressed 


^ _ /2 

V  tanh^  (re/2) 


Sinee  /  is  Oi^nh'f  ,  equation  1(60)  I  ean  be  written  as  the  following  Taylor  series  expansion: 


u)^  /^tanh(re/2)  tanh^  (re/2)  ,  >6 

tanh(re/2)  2^7  8^7  +0[Kh)  . 
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Now,  expanding  the  hyperbolie  functions  on  the  right  hand  side  of  equation|(61)|in  terms 
of  nh  to  Oi^nh'f  ,  and  then  expressing  each  nh  term  as  a  function  of  u  and  h  using  the 
expansion  of  the  finite  depth  dispersion  relation, 


we  obtain 


-\-0[K.h)^ , 


(62) 


6  ^  360  2u;  ^  12 


(63) 


in  agreement  with  equation  I  (5  8)1  Figure  4.2  shows  the  first  (0(l)),  second  (O^nhf), 

and  third  {Oi^nKf)  order  approximations  of  the  cross-shore  wavenumber  compared  to 

the  exact  finite  depth  solution  for  a  normally  incident  wave  with  a  period  of  10s.  Figure 
4.3  shows  the  same  comparisons  for  an  incidence  angle  of  30  degrees. 
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Figure  4.2.  Cross-shore  wavenumber  vs.  depth  for  the  first,  seeond,  and  third  order 
approximations  are  compared  to  the  exact  solution.  Results  are  shown  for  a  10s  period 

wave  at  normal  incidence. 
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Figure  4.3.  Cross-shore  wavenumber  vs.  depth  for  the  first,  seeond,  and  third  order 
approximations  are  compared  to  the  exact  solution.  Results  are  shown  for  a  10s  period 

wave  at  a  30  degree  angle. 


Though  not  a  significant  improvement  for  depths  less  than  20m  where  nh  <l ,  the  third 
order  approximation  does  show  a  notable  improvement  over  the  second  order  for  a  wide 
range  of  intermediate  (40- 100m)  depths.  Note  that  the  third  order  approximation 
eventually  diverges  from  the  exact  result  in  sufficiently  deep  depths;  as  expected  from  a 
Taylor  series  expansion  for  small  nh ,  this  divergence  grows  rapidly  with  depth  compared 
to  the  second  order  result. 

Figure  4.4  shows  the  same  approximations  of  cross-shore  wavenumber  k 
compared  to  the  exact  finite  depth  solution  for  a  10s  wave  in  a  fixed  depth  of  20m,  versus 
incidence  angles. 
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Figure  4.4.  Cross-shore  wavenumber  k  vs.  incidence  angle  6*o .  The  first,  second,  and 
third  order  approximations  are  compared  to  the  exact  linear  solution.  Results  are  shown 

for  a  10s  period  wave  in  20m  depth. 


The  first  order  approximation  does  not  account  for  directional  effects  and  has 
large  errors.  Both  the  second  and  third  order  approximations  accurately  account  for 
directionality.  The  third  order  approximation  is  within  a  few  percent  of  the  exact 
solution  even  for  angles  up  to  60  degrees. 

C.  SHOALING 

The  conservation  of  energy  for  a  wave  train  with  radial  frequency  uo  and 
alongshore  wavenumber  /  may  be  expressed  as 


A 

dx 
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or 


dE  _ 

1  dc^ 

\  1 

1  dn 

dx 

,  dx 

\  / 

K  dx 

(64) 


where  E  =  \A\  is  proportional  to  the  wave  energy,  and  is  the  cross-shore  component 
of  the  group  speed  .  For  arbitrary  finite  depth,  the  (non-dimensional)  linear  dispersion 
relation  is  given  hy  =  k  tanh(KA),  yielding  the  following  analytical  expressions  for 


dx 


dn 

and  — : 
dx 


_  du)  u 

1  nh 

dn  K, 

1 

2  sinh  («;/?) 

(65) 
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1 
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dn  —K,^  dh 

^  sinh^  (inh)  -\-  nh 


(67) 


Substituting  equations  |(65)||(67)  [into  equation  |(64)  [yields 
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In  order  to  compare  the  exact,  linear  energy  equation  with  the  Boussinesq  form 
(•Sd)  we  expand  equation l(ih^ in  a  Taylor  series  for  small  nh  to  ,  which  yields 


dE 

/ 

1 

[rdij  1 

dx 

Ih  2h 

K-l\ 

dx 


(69) 


I  l2  He 

Now,  since  =  |^|  =  AA  (where  denotes  complex  conjugation),  the 

amplitude  evolution  equation  [(5 7) [yields  the  following  energy  equation  (dropping  the  p,q 
indices) 


dE  dA*  ,  dA 

—  =  A - VA  — 

dx  dx  dx 


1 

- h  S 

- 1 - 7 

2h 

1  2  ] 

m 

dx 


(70) 


By  observation,  it  is  immediately  seen  that  equation  69)  and  equation  |(70)  [agree 
to  first  order.  To  compare  the  0(s)  terms,  note  that  (l  +  0(e))  and 


—  (l  +  (9(e)) ,  and  thus  the  0(e)  terms  in  equation  (70)lare  in  agreement  with 

equation  [t69llto  0(e^) . 

The  improved  energy  conservation  properties  of  the  present  model  can  be 
observed  in  Figures  4.5  and  4.6,  which  compare  the  second  and  third  order 


approximations  of  the  energy  flux  coefficient  — 


1  dc 


gx 


^gx  9x 


to  the  exact  linear  finite  depth 


solution.  Again,  results  are  shown  for  a  10s  wave  at  normal  incidence.  Figure  4.5,  and 
for  a  30  degree  incidence  angle.  Figure  4.6. 
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Figure  4.5.  Energy  flux  coefficient 


1  dc 


gx 


dx 


vs.  depth  for  the  seeond  and  third  order 


approximations  are  eompared  to  the  exaet  solution.  Results  are  shown  for  a  10s  period 

wave  at  normal  ineidenee. 
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Figure  4.6.  Energy  flux  eoefflcient 


1  dc 


dx 


vs.  depth  for  the  second  and  third  order 


approximations  are  compared  to  the  exact  solution.  Results  are  shown  for  a  10s  period 

wave  at  a  30  degree  angle. 


Note  that  the  second  order  approximation  diverges  from  the  exact  solution  even  in 
relatively  shallow  water  (<20m),  and  that  this  divergence  is  even  more  extreme  for 
oblique  waves.  The  third  order  approximation  shows  significant  improvement  over  the 
second  order  in  both  the  normal  and  oblique  cases. 

To  illustrate  the  improved  shoaling  characteristics  of  the  third  order  model,  the 
amplitude  evolution  of  a  wave  shoaling  on  a  plane  beach  was  predicted  numerically  using 
a  simple  first-order  finite  difference  approximation  of  the  amplitude  evolution  equation. 
Figure  4.7  compares  the  evolution  of  jA/Aol  (where  Ao  is  the  complex  amplitude  of  the 
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incident  deep-water  wave)  predieted  by  the  seeond  (HB)  and  third  order  model  with  the 
exact  finite  depth  solution. 


AiFitULHle-^fviHud'xi  kr  F  =  A  I  Iti 


Figure  4.7.  Shoaling  of  a  10s  period  wave  normally  ineident  on  a  plane  beach  with  a  1:50 
slope.  Topography  varies  in  the  eross-shore  direetion  only.  Top  panel  compares  the 
seeond  and  third  order  model  results  to  finite  depth  theory.  Bottom  panel  shows  the 

depth  profile. 


The  signifieant  improvement  of  the  third  order  model  over  the  seeond  order  FIB 
model  is  evident  in  the  predicted  amplitudes  in  shallow  water.  For  this  relatively  deep 
(25m)  model  initialization,  the  seeond  order  model  shows  the  expeeted  overshoaling, 
overpredieting  the  wave  amplitude  in  Im  depth  by  28%;  whereas  the  third  order  model  is 
within  1%  of  the  exaet  finite  depth  theory  prediction. 
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V.  SUMMARY  AND  FUTURE  RESEARCH 


In  general,  wave  transformation  on  beaehes  is  well-deseribed  by  one-dimensional 
models.  However,  the  direetionality  of  waves  is  of  erueial  importanee  to  a  variety  of 
nearshore  proeesses  ineluding  infragravity  motions,  alongshore  eurrents,  and  sediment 
transport.  Additionally,  alongshore  varying  bathymetry  eomplieates  nearshore  wave 
evolution.  In  this  thesis,  Boussinesq-type  equations  for  weakly  nonlinear,  weakly 
dispersive  waves  are  derived  for  a  beaeh  with  weak  alongshore  depth  variations. 

We  extend  the  derivation  of  Herbers  and  Burton  (1997)  to  one  higher  order  in 
nonlinearity,  dispersion,  and  bottom  slope  to  derive  a  eoupled  set  of  amplitude  evolution 
equations  for  a  speetrum  of  waves.  The  main  result,  equation  (54),  more  aeeurately 
represents  dispersion,  direetional  spreading,  and  the  strength  of  triad  interaetions  than  the 
HB  model,  and  extends  the  effeets  of  eross-shore  bottom  slope  to  one  higher  order.  New 
effeets  introdueed  at  the  lowest  order  are  alongshore  depth  variations  and  nonlinear 
quartet  interaetions.  Linear  terms  were  verified  by  analytieal  methods,  and  example 
eomputations  in  uniform  depth  and  for  a  plane  beaeh  quantify  some  of  the  improvements 
of  this  model  over  the  lower  order  HB  model.  In  partieular,  it  is  demonstrated  that  this 
higher  order  model  better  represents  the  vertieal  strueture  of  the  flow  field,  notably 
improves  the  approximation  of  the  eross-shore  wave  number  both  in  intermediate  depths 
and  for  larger  direetional  spreading  angles,  and  more  aeeurately  aeeounts  for  energy 
eonservation.  For  a  typieal  swell  period  and  model  initialization  in  relatively  deep 
(eontinental  shelf)  depths,  it  is  shown  that  the  overpredietion  of  shoaling  is  redueed  from 
0(30%)  for  the  HB  model  to  less  than  1%  for  the  present  higher  order  model. 
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Further  work  is  needed  to  verify  other  terms  in  equation  (54)  No  simple 
analytieal  results  exist  for  beaches  with  two-dimensional  depth  variations.  Therefore,  the 
term  that  accounts  for  alongshore  depth  variation  may  have  to  be  verified  numerically 
using  an  appropriate  refraction  model.  Higher  order  cross-shore  bottom  slope  terms  in 
equation  |(54)|  involving  the  square  and  the  curvature  of  the  bottom  slope  are  perhaps 
more  difficult  to  verify.  However,  Chu  &  Mei  (1970)  derived  similar  terms  for  a  Stokes 
wave  in  arbitrary  depth.  Their  expressions  are  rather  complicated  in  terms  of  hyperbolic 

functions.  An  expansion  of  the  Chu  and  Mei  result  for  small  nh  to  Oi^nh'f  should 
render  terms  that  correspond  to  the  cross-shore  bottom  slope  terms  derived  in  this  work. 

Finally,  the  coupling  coefficients  for  triad  and  quartet  interactions  remain  to  be 
verified.  HB  show  that  in  the  limit  of  small  wave  amplitude  and  bottom  slope  the  steady 
solutions  to  their  evolution  equations  exactly  match  the  second  order  bound  waves 
predicted  by  finite  depth  theory  in  the  shallow  water  limit.  It  should  be  possible 
(although  more  algebraically  tedious)  to  extend  this  approach  to  one  higher  order  and 
validate  both  corrections  to  the  triad  coupling  coefficient  and  the  quartet  coupling 
coefficient  in  equation[(54)[ 
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APPENDIX  A 


The  average  horizontal  flow  iu,vj  over  the  water  column  is  given  by 

- J  udz' =  eUi+e^U2+£^u^+o(^e‘^Y 


u  = 


h  +  T]' 


_  -I  _  _ 

v  = -  f  vdz'  =  a(eVi+£^V2+£^Vj+o(£^)). 

h-\-ri'J  '  '  '' 


(71) 


—  h 


The  reciprocal  of  the  water  column  height  can  be  expressed  by  geometric  expansion  as 

+  0(£'). 
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The  integral  term  of  u  is  expanded  (note  u\  is  independent  of  z') 

V  V 
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and  depth-averaged  flow  variables  can  be  defined  to  correspond  with  each  order  of  u 


39 


(72) 
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hJ  h  h 
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Substituting  |(14)  [and  |(22)|into  |(72)|for  U2  and  m,  respectively,  and  integrating  yields 


and 


Uj  U2  ■ 
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^  3 

dx^ 

d^U2 
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’  3 
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dx'dy' 

dh  du^  ^  hrj^  hrj^ 
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(73) 
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Similarly  it  can  be  shown  that  the  depth-averaged  alongshore  flow  variables  are  given  by 


(75) 


3  dxdy ’ 
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and 
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APPENDIX  B 


10/ 

cj)^ ;  M,.  =  — - ,  V,.  =  — -  .  To 
dx  dy  \ 

make  necessary  substitutions  to  the  depth-averaged  governing  equations,  these  functions 
must  be  expressed  in  a  depth-averaged  form.  For  convenience,  define  to  be  the 
average  value  of  0,  from  the  bottom  to  the  still  sea  surface 


-  1 


(78) 


—  h 


Since  u\  and  vi  are  independent  of  z,  so  is  (j)^.  Applying  Leibniz’s  Integration  Theorem, 


■  _  1  r  d4>2  j  _  d(j)2  ^  (^2  ^2  (  ^))  dh 

h^.  dx 

-h 


dx 


dx 


(79) 


To  eliminate  the  explicit  boundary  value,  02  (~^)  >  from  equation  [79u  recall 


Bernoulli’s  equation 


d(b  M  -|-v  -|-iv 

p  = - z. 

dt  2 


(80) 


Substituting  the  perturbation  expansions  (4b-e)  and  the  definitions  of  0 ,  and  collecting 
terms  of  like  order,  we  obtain 
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Hence,  using  equation  [(7  8)[ 


d(j)2  _  d^Uy  uf 

dt  3  dxdt  1 


and 


so  that 


502  {—h)  _  5^Mj  wf 

dt  2  dxdt  2 


^(02-02H)) 


d\ 
6  dxdt 


Integrating  in  time,  setting  the  arbitrary  integration  constant  to  zero. 


02  ^2  ( 


duy 
6  dx 


Thus,  equation  [79)  [simplifies  to 


—  50,  h  du,  dh 

M,  = - ^  — £ - 

dx  6  dx  dx 


Similarly  it  can  be  shown  that 


502  h  du^  dh 
dy  6  dx  dy 


Applying  the  same  procedure  to  the  next  higher  order  yields 


—  50,  nh  d^u, 
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dy  3  dxdy  ^ 
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(83) 


(84) 
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APPENDIX  C 


Substitution  of  the  assumed  waveforms  for  0j  (equation I (44)b  and  (equation 
1(45))  and  the  solution  for  772  {(pi  —  o)  derived  in  Appendix  D  (equation  l(95)|  into  the 


governing  equation  [(43 )| yields  the  following  expressions  for  eaeh  of  the  terms: 
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APPENDIX  D 


To  determine  eonsider  the  surfaee  boundary  eondition 
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Substituting  the  perturbation  expansions  (4a-d)  and  scaling  relations  |(5)|  into  equation 
1(92)  and  collecting  terms  of  like  order  yields 
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Now,  substitute  the  assumed  lowest  order  wavefield  |(44)|and  |(45)|  (dropping  the  primes 
to  simplify  the  notation)  to  express  each  term  involving  r]„  u„  and  v,  in  terms  of  the 
spectral  amplitudes  : 
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Finally,  integrating  equation  |(94)|  with  respect  to  time  and  setting  the  arbitrary  integration 
constant  to  zero. 
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